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ABSTRACT 

We report two different types of backflow from jets by performing 2D special relativistic hydrody- 
namieal simulations. One is anti-parallel and quasi-straight to the main jet (quasi-straight backflow), 
and the other is bent path of the backflow (bent backflow). We find that the former appears when 
the head advance speed is comparable to or higher than the local sound speed at the hotspot while 
the latter appears when the head advance speed is slower than the sound speed at the hotspot. Bent 
backflow collides with the unshocked jet and laterally squeezes the jet. At the same time, a pair 
of new oblique shocks are formed at the tip of the jet and new bent fast backflows are generated 
via these oblique shocks. The hysteresis of backflow collisions is thus imprinted in the jet as a node 
and anti-node structure. This process also promotes broadening of the jet cross sectional area and 
it also causes a decrease in the head advance velocity. This hydrodynamic process may be tested by 
observations of compact young jets. 

Subject headings: hydrodynamics - galaxies: jets - methods: numerical - relativity 



1. INTRODUCTION 

Recent observations of jets in active galactic nuclei 
(AGN) clearly tell us their morphological charac- 
ters. On VLBI scales, the MOJAVE (Monitoring 
Of Jets in Active galactic nuclei with VLBA Ex- 
periments) program clearly shows us a number of 
actual internal inhomogencous structures, and a 
growth of the jet's cross-sectional area along its axis 
jhttp: //www. physics. pu rdue . edu/ ~ mlist er / M O J AVE / index . 
iLister et al.l (2009)). Their cross sections are not con- 
stant nor internally uniform. The importance of jet's 
internal structures has also been used to explain multi 
frequency spectra (jGhisellini et al.l l2005h . From the 
standpoint of hydrodynamics, however, the generation 
of internal structures in jets and what determines the 
cross sectional area of jets during their evolutions have 
not yet been fully investigated. 

Most prior work has focused on Kclvin-Hclmholtz in- 
stability (KHI) caused by the velo city shear between 
a jet a nd surrounding matt er (e.g., lHardee fc Normanl 
V iPerucho et al. |2006f l). The instability is known 



gas, respectively. Although the existence of backflows is 
well known to previous studies by numerica l simulations 
of th e jet propagation ([Norman et al.lll98 2; Mart i et al 



19971: iRosen et al.1 119991 : lAlov et all 120001 : iZhang et ~ 



2003( 1. the importance of the backflow in terms of dy- 



as a trigger of diamond-shape internal structures caused 
by oblique shocks inside the unshocked jet (e.g, Nor- 
man et al. 1984). Once the unshocked jet is per- 
turbed by the KHI in a cocoon which surrounds the jet 
(jBegelman et al.l 119841 ) . then an oblique shock appears 
inside the unshocked jet and the velocity component per- 
pendicular to the jet axis is canceled out thr ough this 
oblique shock (|Komissarov fc Falld H997. 199J). 

In this Letter, we newly show the importance of back- 
flow. A backflow is usually generated at the hotspot of 
the light jet, where a light jet is a jet with mass den- 
sity (pj) lower than that of the ambient gas (p&), where 
the subscripts 'j' and 'a' stand for the jet and ambient 
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namics is not well recognized. The backflow dynamics 
are usually ignored since the backflow is subsonic flow. 
For the first time, it is shown that backflow is essential 
jhn^prming the internal structure and the growth rate of 
the cross-sectional area of jets. 

2. NUMERICAL SET-UP 

We utilize the axisymmetric two dimensional special 
relativistic hydrodynamic code in cylindrical coordinate 
(r x z) developed by one o f the a uthors (AM: see for 
details, iMizuta et all (|2004l . 12006( 1). Using the PPM 
method and TVD Runge-Kutta time integration, we 
achieve 3rd order accuracy in terms of space and time. 
Since we do not include any radiative process, the equa- 
tions solved in this study are scale free. The ideal gas 
equation of state (p = (7 — l)pe) is adopted, where p is 
pressure, e is specific internal energy, and 7 is specific 
heat ratio which is fixed to 5/3. The local sound speed 

is derived as c s = y/ jp/ ph = \J 67(7 — 1)/(1 + ej). 

The boundary condition at the cylindrical axis (r = 
0) is the reflective condition. The reflective boundary 
condition is assumed at z = except for the first 10 
grid points, i.e., < r < lRb where a continuous jet is 
imposed. The Rb is the radius of the injected jet from 
the computational boundary (z = 0). A zero gradient 
free boundary condition is assumed at z = 100 Rb and 
r = '"max- The unit of the space scale is the beam radius 
of the injected jet (Rb). The unit of the time scale is 
Rb/c, where c is speed of light. 

In this study, five cases with different jet velocities (vj) 
have been examined these are Vj = 0.9c (model J09), 
vj = 0.8c (J08), v 3 = 0.7c (J07), Vj = 0.6c (J06), and 
Vj = 0.5c (J05). Corresponding Lorentz factors (r) are 
2.3 (J09), 1.7 (J08), 1.4 (J07), 1.3 (J06) and 1.2 (J05), 
respectively The jet is constantly injected from one of 
the computational boundaries. Other model parameters 
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are fixed as follows. The density ratio is fixed to be 
r\ = Pj/pa = 0.1. The initial mass density profile of 
the ambient matter p a is assumed to be uniform. Mass 
densities are normalized by the ambient gas, i.e., p a = 1. 
Both the jet and the ambient gas are assumed to be cold, 
i.e., £j/c 2 = 1CP 2 and e a /c 2 = 10~ 6 , where e is specific 
internal energy Thus specific enthalpy (h = 1 + e/c 2 + 
p/p) is ~ 1. Corresponding sound speeds of the jet and 
ambient gas are c s j ~ 0.16c and c s . a ~ 1.6 x 10 _3 c, 
respectively. The velocity vector of the injected jet is 
parallel to the cylindrical axis. We intend to do extensive 
studies with wider parameter space in the near future. 

The numerical runs are performed until the head of 
the jet reaches the boundaries at z = 100i?6. Since 10 
uniform grid points per a Rb are given in all computa- 
tional domains for both radial and z directions, there 
are in total 1000 grid points in the z direction. Since the 
aspect ratio of the distances between the propagation di- 
rection and lateral direction is different in each model, 
we set different computational domains (r max ) in r by 
the models so that the sideways expanding shock is in- 
side the computational domain in the last stage of the 
computation. 

3. NUMERICAL RESULTS 

Here we mainly discuss the results of J09 and J05 in 
order to study differences in backflow dynamics more 
clearly. The other three cases will be briefly summarized 
later on. Figure [T] shows the mass density contours at 
the final stage of each model (model J05: t = 1845 Rb/c, 
model J09: t = 255 Rb/c). Significant differences are 
found between the two models. They are characterized 
as follows. 

(PI) Significant nodes/anti-nodes structure is generated 
along the jet in J05. The cross sectional radius 
varies along the jet in J05, while the cross sectional 
radius is almost constant along the jet in J09. 

(P2) The radii of the jet at the terminal shock are about 
1.5 Rb in J09, whilst the radii is 5 Rb in J05 which 
is wider than that in J09. 

(P3) The aspect ratio of the surrounding bow shock 
appearing in J09 is narrower than the one in J05. 

In both models, we observe the backflow which begins the 
hotspot at the head of the jet. The temperature in the 
hotspot is quite high, since all the jet's kinetic energy is 
converted to thermal energy through the terminal shock. 
The local sound speeds in the hotspot are c s j ls ~ 0.5c 
(model J09) and c Sj hs ~ 0.3c (model J05). In order to 
clarify the physical reasons for these differences in the 
density map, below we examine the time evolution of 
each model in detail. In advance, we stress that the ve- 
locity field maps are the key to understanding the above 
difference because backflows play an important role in 
these maps. 

Figure EJa) shows the contour of the absolute velocity 
(color) and velocity vectors (arrows) of the model J09 
in its early phase (at t = 30 Rb/c). The jet (red re- 
gion in contour in Fig|Sfa)) ends at the terminal shock 
(at z = 13 Rb), forming a hotspot. We find that the 
jet is surrounded by a layer which is made of shocked 
matter escaped from the hotspot. The velocity of this 



flow is quite small and it sometimes shows a negative 
value measured from the observer frame, i.e., v z < 0. 
This is because c Sj h s and the head propagation speed 
Vh ~ (14i?&)/(30i?b/c) ~ 0.5c are very similar, where 
Vh is the speed of the contact discontinuity between the 
shocked jet and the shocked ambient gas along the jet 
axis. The path of this narrow flow is also almost straight 
near the head (quasi-straight backflow) in the late phase 
of evolution, (see Fig. [2fb)). Since the straight backflow 
is unstable for the KHI between the jet and the backflow, 
some vortices are generated along the jet. 

Figures EJc)-(f) show the evolution of the absolute ve- 
locity contours for model J05. It is found that the path of 
the backflow of model J05 differs dramatically from that 
of model J09. At t = 105R b /c, v h ~ (13R b )/(105R b /c) ~ 
0.1c. The path of the backflow from the hotspot is not 
straight but bent (bent backflow) near the head, (see 
Fig. EJd)). The speed of the bent backflow near the 
hotspot is ~ 0.25c. This bent backflow compresses and 
piles up the shocked ambient gas. Then, the backflow 
makes a half loop path. Finally, the bent backflow ver- 
tically collides with the main jet from the side, (see Fig. 
0(d) at z = 7 Rb). 

As shown above, the transition of backflow dynamics 
is caused by different head advance velocity. The head 
advance velocity decreases from model J09 to model J05. 
The intermediate models (J08, J07) show quasi-straight 
backflow paths, while the backflow of the model J06 ex- 
hibits bent path when each head reaches z w 12R b , see 
Fig [31 

4. PHYSICAL INTERPRETATIONS 

In this section, we explain the physical reasons for the 
significant difference between the two models. According 
to our simulations, the backflow dynamics are the most 
important feature in understanding these differences and 
we have found mainly two features. One is the difference 
in the ratio of Vh to c Sj hs- As we have shown, is com- 
parable to c s .hs in model J09, whereas the head propa- 
gation speed is slower than local sound speed in model 
J05. Since the backflow starts the hotspot, the speed of 
the backflow at the frame of the hotspot is close to c s .h s - 
If Vh is slower than the local sound speed, the backflow 
expands sideways. In the other case, i.e. model J09, the 
backflow goes directly backward and flows anti-parallel 
flow to the main jet, see a schematic figure of the two 
types of the backflow path (top panel of Fig. [4|). By the 
comparison of Fig.[2^a) and Fig.[2]^d), it is found that the 
jets propagated at the same distance from the two mod- 
els but along different paths of the backflow, i.e., quasi- 
straight backflow and bent backflow. The bent backflow 
is generated when the condition c s .h s /^h < 1 — 2 holds. 
For Cg.hsA'h ~ 2, the backflow is clearly separated from 
the jet. 

The other reason is the emergence of the oblique shock 
at the tip of the jet-terminal which is shown in Fig. Q] 
(the red line in the top panel). The bottom panel of 
Fig |4] shows log-scaled density (top) and pressure (bot- 
tom) contours of the same model and at the same time 
of Fig. M (e), i.e., model J05 at t = 270R b /c. The 
density and pressure discontinuity can be clearly seen 
at 15 < z/R b < 20, \r\ ~ lR b . This is an oblique 
shock and is indicated by the d ashed line in the c on- 
tours (Fig. gl), see also Fig. 8(b) ijMizuta et al.ll2004fl in 
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which the oblique shock and fast backflow emerging thor- 
ough the shock can be clearly seen. When the oblique 
shock appears, the backflow comes out not only from the 
hotspot but also from the oblique shock, as shown in the 
shematic figure in Fig[4j A weakly dissipated flow passes 
through such an oblique shock and it becomes bent back- 
flow. The speed of the bent backflow is ~ 0.5c. Although 
little attention has been paid to such a weakly dissi- 
pated flow, it is essential for producing backflows. This 
wea kly dissipated flow was seen in some prev ious studies 
i.e.. lSaxton et alJ (|2002f ); lMizuta et all (|2004f ) . Since this 
flow does not pass through hotspots, the flow directly be- 
comes a backflow, keeping its high velocity. This flow is 
separated from the main jet and draws a half loop in the 
cocoon. As a result, a new oblique shock emerges near 
the terminal, when the bent backflow vertically hits the 
unshocked main jet. 

PI and P2 are understood as follows. Above, we eluci- 
date that the ratio "Uh/c s ,hs divides the types of backflows. 
Below, we explain how PI and P2 are achieved. When 
the bent backflow laterally collides with the unshocked 
jet, the jet is squeezed by the dynamical pinch. Thus 
oblique shocks are generated inside the jet. Just ahead of 
the squeezed node, the jet starts to expand like a nozzle. 
This process for generating oblique shocks is completely 
different from the previously known one in which internal 
shocks are caused by a pressure fluctuation in the cocoon 
(e.g., Norman et al. 1982). Since the terminal shock is 
close to the collisional point, this oblique shock does not 
work very well to rcconfinc the jet. On the contrary, the 
cross section of the jet increases at its terminal as shown 
in the schematic figure (Fig. 3]) . The wider the cross sec- 
tion becomes, the slower the jet advance speed becomes. 
Furthermore, a pair of oblique shocks represented by the 
red line reaches the terminal of the jet and connects with 
the reverse shock. Once the bent backflow appears from 
the hotspot, it is easy to generate these oblique shocks, 
and new bent backflows intermittently emerge from them 
as discussed above. The resultant structure therefore has 
significant node/anti-node structures which are the hys- 
teresis of the backflow attacks as pointed out in PI and 
P2. 

In J09, Vh is almost constant throughout the whole cal- 
culation. In J05 Vh is also constant up to t ~ 200 Rb/c. 
After that, the head propagation velocity be gins to decel- 
erate a s is shown in the numeri cal studies of Scheck et alJ 
(|2002j ) and lMizuta et all (|2004l ). A large vortex is formed 
near the hotspot in the simulations. The large vor- 
tex causes the decelerat ion of the head propagation. 
iKawakatu fc Kinol (|2006f ) analytically obtain taking 
the growth of head cross sectional area into account and 
it reproduces these numerical results well. Thus, P3 can 
be simply understood in terms of previous studies, i.e., 
the deceleration of v^. Finally, it should be noted that 
the bent backflow is indirectly linked with the generation 
for a large vortex, since a weak dissipated flow through 
an oblique shock is a key for the appearance of the vortex 
(jMizuta et al.ll200l . 

So far, we fix 7 = 5/3 for simplicity. This is inappro- 
priate when the specific internal energy is in a relativis- 
tic regime (e/c 2 > 1), and it could be achieved in the 
hotspot. Since the specific heat ratio should be smaller 
than 5/3 in the relativistic regime, the actual backflow 
velocity from the hotspot is slightly slower than our es- 



timation. Therefore, the transition velocity of backflow 
dynamics is also smaller. In order to check the validity 
of our discussions thus far in the real system, we exam- 
ine the case of fixed 7 = 4/3 which corresponds to the 
relativistic limit. Three injection velocities, vj = 0.4c 
(model J04-43), vj = 0.5c (model J05-43), and vj = 0.9c 
(model J09-43) are studied. Figure[5](a)-(c) show the ab- 
solute velocity of these models in the early phase of the 
evolution. These simulations also show that there are 
two types of backflows, the quasi-straight path (models 
J05-43 and J09-43) and bent path (models J04-43). It 
is verified that the transition velocities are Vj ~ 0.4c, 
whereas Vj ~ 0.6c in the fixed 7 = 5/3 case. That is, the 
true transition would be a certain value between them. 
Therefore we conclude that our simplification does not 
change the essence of our new findings. 

5. SUMMARY AND DISCUSSION 

In this work, we show two types of backflow dynamics, 
quasi-straight path and bent path. The difference can be 
clearly seen in the early phase of the dynamics, (see Figs. 
[2] (a) and (c)). When the head propagation velocity is 
smaller than local sound speed at the hotspot, a bent 
backflow appears. The bent backflow goes backward and 
sometimes attacks the jet from the side. When the jet 
is dynamically attacked by the bent backflow, the cross 
section of the jet increases, making an oblique shock. At 
the terminal of the jet, an oblique shock and the ter- 
minal reverse shock connects, making a weak dissipated 
flow into the cocoon through the oblique shock. Such 
backflow also becomes a bent one and appears repeat- 
edly. The transition of two types of backflow dynamics 
occurs by decreasing the head propagation velocity. 

Compact Symmetric Object (CSO) would be one of 
the best targets for testing our findings, since CSOs are 
newly born AGN jets with radio lobes. For instance, 
one of the smallest CSO radio lobes 3C 84 lAsada et al.l 
(|2006D with the size of ~ 10 pc may be a good target to 
explore the backflow physics investigated in this work. 
Observations with unprecedented high spatial resolution 
which will be provided by VLBI Space Observatory Pro- 
gramme 2 (VSOP-2) in the near future appear to provide 
a promising means of exploring these backflows. 

We add brief comments on the relevance to gamma-ray 
bursts (GRBs). Hot jets are expected in GRBs and it is 
intriguing to see how backflows act when jets are hot. 
Since the GRB jet from collapsars is formed deep inside 
of the progenitor, the jet should drill the high density 
stellar envelopes to reach the progenitor surface (e.g., 
ILazzati et alJ (|2009| )). For example, the bent path of the 
backflow from the head of the GRB j et is indeed seen 
in Figures 2, 3 and 4 in I Zhang et al.l (|2004D , during its 
propagation in the high density progenitor envelope. 
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Fig. 1. — Density contours of the models J05 (left) and J09 (right) at the end of simulations. In both cases internal structures can be 
seen inside the jet. The jet radius increases at anti-node and decreases at the node for the case J05, while the radius of the jet is almost 
constant along the jet for the case J09. 
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Fig. 2.— Absolute velocity contours of the model J09 at (a) t = 30 and (b) 255 R b , and of the model J05 at (c) t = 30, (d) 105, (e) 270, 
and (f) 1845 Rb- The arrows show the velocity vectors. The color bars are common for (a) and (b), and for (c)-(f). The backflow is almost 
straight for the model J09. Whereas the backflow is bent for the model J05. A bent backflow attacks the jet from the side (d)-(f). 
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Fig. 3.— Same as Fig. E] but for models (a)J06 (t = 75 R b /c), (b)J07 {t = 52.5 R b /c), and (c)J08 (t = 37.5 R b /c). The jet propagates 
about 12R b in each models, see also Fig. 2(a) and (d) for the comparison. The path of the backflow is quasi-straight for models J07 and 
J08. On the contrary the path of the backflow is bent for model J06. 



8 



Mizuta et al. 



Vu > c 



s'hs 



hotspot 



contact discontinuity " 

backflow 



bow 



jet ^ 




cocoon backflo^ s ~*' 

shock ed ambient gas / g^^ e 



v h <c 



J05 t=270Rb/c 




- jet 



-5 - 



10 



oblique shock 



10 



log(p) 

-0.5 
-1.5 
-.2, 

log (p) 



- 



-3 
-4 
-5 
-6 



Fig. 4. — A schematic figure of the path of the jet and backflow (top). Two cases are shown. The head speed is faster than the local 
sound speed at the hotspot (upper top) and the head speed is slower than the local sound speed at the hotspot (lower half). The path 
of the backflow is quasi-straight for the former case. The path of the backflow is bent for the latter case. The backflow from the hotspot 
thermally expands sideways, and becomes bent backflow, interacting with the shocked ambient gas. An oblique shock sometimes appears 
(red line) at the terminal of the jet which is another way to trigger bent backflow. The bottom panel shows log-scaled density and pressure 
contours of the same model and at the same time of Fig. [21(e), i.e., model J05 at t = 270-R;,. Near the terminal shock which is at z = 20-R;, 
and perpendicular to the z axis, a discontinuity indicated by the dashed line can be seen. This corresponds to an oblique shock shown by 
the red line in the top panel. The backflow mainly comes through this oblique shock at this snapshot, see also Fig. |2](e) in which a fast 
backflow starts through this oblique shock and it becomes a bent backflow. 
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Fig. 5.— Same as Fig. |2] but for models (a) J04-43 (t = 127.5 R b /c), (b) J05-43 (t = 82.5 R b /c) , and (c) J09-43 (t = 30 fl 6 /c). These 
simulations are done with fixed specific ratio 7 = 4/3. The path of the backflow is quasi-straight one (model J05-43 and J09-43). On the 
contrary the path of the backflow is bent one for model J04-43. The transition of backflow property occurs at slower injection velocity than 
that for the models done with fixed specific heat ratio 7 = 5/3. 



